Operator fidelity susceptibility, decoherence and quantum criticality 
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The extension of the notion of quantum fidelity from the state-space level to the operator one 
can be used to study environment-induced decoherence. State-dependent operator fidelity sucepti- 
bility (OFS), the leading order term for slightly different operator parameters, is shown to have a 
nontrivial behavior when the environment is at critical points. Two different contributions to OFS 
are identified which have distinct physical origins and temporal dependence. Exact results for the 
finite-temperature decoherence caused by a bath described by the Ising model in transverse field are 
obtained. 



I. INTRODUCTION 



Generically the interaction between a quantum system 
and its environment results in decoherence, and may lead 
the system to experience the so-called quantum-classical 
transition Jj • For this reason the decoherence process is 
regarded as the main obstacle for jthe implementation of 
quantum information processing [3] . Generally speaking 
the properties of the environment may strongly affect its 
decohering capabilities [3, Ij, [a, 0, @1 ■ This implies that 
a quantum system can be regarded as a probe to extract 
useful information about the coupled environment, e.g., 
the quantum phase transitions (QPTs). Very recently 
a nuclear magnetic resonance (NMR) experiment of this 
kind has been performed [^. There a single qubit was 
used as a probe to detect the (precursors of the) quantum 
critical point of the coupled environment. 

The relationship between the decoherence of a central 
spin and the QPTs of the coupled environment can be 
established through the notion of Loschmidt echo (LE) 
Q. It is well known that LE can be exploited to mea- 
sure the stability of a quantum system to perturbations 
0) [l3, \lM- ^^ ^^f- 0' ^ central spin coupled with an 
Ising spin chain in a transverse field was considered; the 
authors found that the decay of LE is enhanced by quan- 
tum criticality. The connection between LE and quan- 
tum criticality has been further clarified in Ref. [0] . There 
the authors showed that the LE enhancement holds true 
just for short decay times (gaussian regime). A genuine 
signature of criticality, on the other hand, can be recov- 
ered by studying the asymptotic (large time) behaviour 
of LE as function of the environment coupling parame- 
ter Q. More recently, the averaged LE over all states 
on Hilbert space with a Haar measure, called operator 
fidelity, was proposed to study quantum criticality [12| . 
This proposal is a direct extension to the operator level 
of state-space quantum fidelity. This latter notion re- 
cently attracted a lot of attention as a new tool to ana- 
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lyze quantum criticalit y, b oth a zero-temperature and at 
finite temperature ^,M, M,M,\l^^- 

Whereas in Q the initial state of the spin chain was 
assumed to be the ground state, in this paper we aim 
at considering a more general situation: the initial state 
of the environment is a mixture of the eigenstate of its 
Hamiltonian. The Gibbs thermal state is special instance 
of this setup. With this assumption the study of decoher- 
ence naturally leads to a state-dependent generalizations 
of the operator fidelity suggested in [I4I . 

In this paper we will discuss the operator extensions of 
fidelity on general grounds. In order to do that one has 
simply to notice that finite-time quantum evolutions cor- 
respond to unitary operators that themselves belong to 
an Hilbert space i.e., the Hilbert-Schmidt one. It follows 
that some of the results obtained in the fidelity approach 
for quantum states lifts immediately to the operator level 
we are now interested in. We will show that the corre- 
sponding operator fidelity susceptibility (OFS) contains 
two different contributions. These two terms arise from 
variation upon parameter change of energy levels and 
eigenstates respectively and have very different tempo- 
ral behaviour. We will then get a general expression for 
OFS for models with a factorization structure typical of 
quasi-free models. Finally, we will exploit OFS for study- 
ing an environment described by an Ising chain with a 
transverse field that can be driven to quantum criticality. 

This paper is organized as follows: In Sec. |TT1 we give 
a general description of the models and the relations be- 
tween decoherence of the central system and the opera- 
tor fidelity of two time evolutions of the environment. In 
Sec. mil we give generalities about operator metric and 
fidelities and analyze the general OFS. In Sec. IIV Al we 
consider the models with a factorization structure, and 
give a general expression for the OFS, and then study 
the Ising model with a transverse field, get the exact so- 
lution for the OFS. Sec. |V] contains the conclusions and 
outlook. 



II. DECOHERENCE AND OPERATOR 
FIDELITY 



III. METRICS OVER MANIFOLDS OF 
UNITARIES 



As mentioned in the introduction, an important phys- 
ical motivation for the operator fidehties we are going 
to analyze in this paper is provided by the decoherence 
process. To see this clearly, let us consider a purely de- 
phasing coupling of a system with its environment. The 
interaction Hamiltonian has the form 



Hj 



^|n)(n| 



Br, 



(1) 



The entire Hamiltonian is H = Hq + Hj , where Hq con- 
sists of system Hamiltonian Hs — X^n E!n\n){n\ and bath 
Hamiltonian Hb- The Hg eigenstates {|?^)}, play here 
the role of preferred pointer states [l[ . 

With the initial state p{0) = |i/'(0))(V'(0)| (g) pb, the 
reduced density matrix of the system at time t can be 
written as 



PS 



(t) = Yl "^""^ 



,-iiE„-E,^)t 



n){m\TTB[pBViit)Vnit)], 



(2) 
where c„ = (n|V'(0)), K(0 = exp [-it{HB + B„)] 

Obviously, in the generic case ps evolves from ini- 
tial pure states to mixed states. The decay of the off- 
diagonal elements of ps means a reduction from a pure 
state to a classical mixture of the preferred pointer states 
— quantum-classical transitions. The temporal behav- 
ior of the off-diagonal element is decided by two factors, 
one is CnC^, just relating to initial state, and the other 
T^'^B[pBV^{t)Vn{t)] is unrelated with the initial state of 
the system, reflecting the dephasing effect induced by en- 
vironment. The latter can be considered as a fidelity for 
two operators Vm and Vn- which is defined by 



Fp{V^,Vn)=\{Vrn,Vn)p 



(3) 



where the inner product {Vm,Vn)p = TY[pV^Vn]. This 
fidelity can be obtained from the density matrix of the 
central system and encodes information about the bath 
state. Notice that one might choose B„ = A„i?; in this 
case different Vn simply correspond to different values of 
the coupling strength A in front of the "perturbation" B. 
This is the scenario we mostly have in mind in this paper 
(see Sec. HVl) . 

Let us now specialize to a 2-level system coupled with 
a bath. If interaction is weak enough, the two effective 
Hamiltonians are slightly different. For an instance in 
Sec. IIVI we will consider the important case where the 
initial state of the bath is a thermal equilibrium state 
PB = exp{-PH)/Z{l3) [Z{(3) = Tr[exp(-/3iJ)] is the par- 
tition function, T = /3^^ is the temperature of the bath]. 
In such a case | Fp^ (14, V2) P will become the Loschmidt 
echo for /3~^ = (for non degenerate ground stat^and 
it coincides with the operator fidelity defined in [12] for 
(3 = 0. 



In this section we discuss operator fidelity from a rather 
general mathematical point view. Let 7i be a quantum 
state-space and a {U\} C U{H) a smooth family of uni- 
taries over Ti. parametrized by elements of a manifold M . 
Given a state p e S{n) = {p S C{n) ] p > 0, trp = 1} 
one can define the following hermitean product over C{7i) 

{X,Y)p:=TTipX^Y) (4) 

If p > then ^ is non-degenerate and \\X\\p := 

y^{X,X)p defines a norm over C{H). In general, if kerp 7^ 

{0}, ]j»]|p is just a semi-norm (if the range oiX is included 

<ll^ll 
= 1. 

Definition The p-fidelity of the operators X and Y is 

given by 



in (suppp)-*^ then \\X\\p = 0) Notice that i) \\X\\p 
and ii) unitaries U arc normalized i.e., {U, U) p 



Tp{XX)--=\{X.Y)^ 



(5) 



It is immediate yet important to realize that for non 
full-rank states p, having two unitaries with p-fidelity 
one does not imply their equality (up to a phase). 
Indeed, from the Cauchy-Schwarz inequality on has 
!Fp{X,Y) < \\X\\p\\Y\\p, in particular for unitaries U 
and W, J'p{U,W) < 1. One has that J'p{U,W) = 1 ^ 
C^jsuppp = W^jsuppp where suppp := (kerp)-'-. 

In order to unveil the operational meaning of the above 
definition it is useful to recall the following fact. If ]^p) e 
H (^H is a, purification of p = J2iPi\^){^\ i-^-, l^p) — 
Ei VP'i\i) <» K> then {X, Y)p - (Vpl(Xt ® l){Y^l)\^p). 
The operator scalar product ^ can be seen as a scalar 
product of suitable quantum states of a bigger system. 
This simple remark shows that operator fidelity ([5]) quan- 
tifies the degree of statistical distinguishability between 
the two states 

\^iA)):=iA®t)\^p)iA = X,Y). 

When U and W denote unitary transformations the 
quantity ([5]) has an interpretation as visibility strength in 
properly designed interferometric experiments [25| . Yet 
another kind of operational relevance of the operator fi- 
delities [5] in the context of decoherence will be discussed 
in HI 

Now, following the differential-geometric spirit of [20| . 
we are going to consider the operator fidelity ([5]) between 
infinitesimally different unitaries. The leading term in 
the expansion of ([5|) will define a quadratic form over the 
tangent space of the manifold U{TL). For full rank p that 
quadratic form is a metric. The following proposition 
shows that and its proof it is just a direct calculation 
analogous to the one performed at the state-space level 

0. 

Proposition Let {U\] C U{TL) be a smooth family of 
unitaries over Ti. parametrized by elements A of a mani- 
fold M. One finds 



Tp{UxXx+sx)^l 



5X^ 



XpW 



(6) 



where, if U' — dU/dX one has 



(7) 



The quantity above wiU be referred to as operator fi- 
dehty susceptibihty (OFS). In the following we will be 
discussing three different situations fitting in the frame- 
work we have set. The first example (Ex. 0) will show 
in which precise sense the operator approach here un- 
der investigation includes the ground-state fidelity one. 
In Ex. 1 we notice that when the family of unitaries is 
a one-parameter group generated by an Hamiltonian H 
then the metric is nothing but the variance of H. Fi- 
nally, in Ex. 2 the most important case of a fixed-time 
family of evolutions generated by a parametric family of 
Hamiltonians is analyzed. 

Example Let {H\} a family of isodegenerate quan- 
tum Hamiltonians with spectral resolutions 

H{X) = Y,EnW\M^)){M^)\- 

n 

The adiabatic intertwiners are unitaries such that 
0(A,Ao)|V^„(Ao)> = |^„(A))Vn. Let p = |V^o(Ao)>(^o(Ao)| 
(lV'o(Ao)>=Ground State of if (Aq)) and U{X) = 0(A, Aq) 
then 

{U{X),U{X'))p = (V'o(Ao)|Ot(A,Ao)0(A',Ao)|V'o(Ao)) = 
(V'o(A),^o(A')>. (8) 

This shows that the ground state fidelity is a particular 
instance of the setting we are discussing. 



Example 1 Let {Ut 



-itH 



}t and p 



Then 



U' = -iHU{t) whence Xp{t) = (0|if^0|) - \{(f>\H[ 



Example 2 {Ux = 



•,-itH^ 



}\ and p — t/d {d := 



dimTi). From time-dependent perturbation theory one 
obtains 



U' ^-iUx f dTe""^H 
Jo 



irHx fT'^-irHx 



(9) 



Introducing the superoperator Lh :~ [H,'] over 
C{Ti) one sees that, if |n)'s and the E^s denote the 
eigenvectors and eigenvalues of H{X), En — E,n and 
^n,m '■= vd\n){m\ are the eigenvalues and (normal- 
ized) eigenoperators of Lh, respectively. One has 
{U',U') = {J^ dre'^^" {H')J^ dae""^^ {H')). If P = 

J2n \^n,n){^run\ IS the projection over the kernel of Lh 
and Q = 1 - P one finds e*"^^" 
the above expression becomes 



Qe'^-^«Q + P, therefore 



(c/',f/'> = (Q' 



^itLh _ \ 



^H 



Q{H'),Q- 



airLH 



^H 



'-QiH')) 



t^{P{H'),P{H')) ^\\FtQ{H')f +t^P{H')f (10) 

(notice that Lh is invertible, by definition, on the range 
of Q) where Ft{x) := s\n{xt/2)/{x/2). 

The first term in the last line of ([T(])) can be rewrit- 
ten as {H'\QF^{Lh)Q\H'). Similarly the term {U,U') 



in the metric can be written as {1, L dre^'^^" {H')) — 
J*dT(e-''^"{l),H') = t{l,H') = d'^TrH'. Putting all 
together and expanding over the basis |^n,n) of Lh- 
eigenoperators one finds the explicit result 



Xp(A) = :, > I 



1 x-^ .(nlH'lm) 



|2t72 



d ^ En — E„ 



Ft{En~E.rn) (11) 



t\W\{n\H'\n)\^ ~ ^TrHr) (12 



The key property here is that 



lim t'^F^{x) = 2t:6{x). 

t — ^oo 



This asymptotic delta function is responsible for the large 
contribution to p^ given by small £"„ — Em This shows 
that all (quasi) level crossings in the spectrum of H might 
lead to a analyticity breakdown in Xp- This has to be 
contrasted with the GS fidelity studied e.g., in [l^ where 
just (quasi) level crossings in the GS play a role. In this 
sense it is clear why the OFS p^ is a more powerful tool 
than the corresponding state-space analogue. 

An important generalization of (fT2)) can be obtained 
by considering p commuting with H. In this case one 
obtains an expression analogous to (fT^ with the diagonal 
elements p„^„ = {n\p\n) suitably inserted (see Sec. lIII Al) . 



A. Splitting Xp 

Now we derive an alternative form for the OFS ([7]) for 
the case of Ex. 2 discussed above. This form will make 
even more manifest the different contributions to opera- 
tor suceptibility arising from eigenvalue and eigenvector 
variation with the control parameter A. An unitary op- 
erator can be written in diagonal form 



C/(A) = ^«.(A)P.(A). 



(13) 



where the m^'s are the eigenvalues of U satisfying \ui\ = 1 
and Pi — \ipi){il)i\ is an one-dimensional projective oper- 
ator. Both eigenvalues and eigenstates are parameter 
dependent. What we would like to do first is to sepa- 
rate the contributions of these two kinds of parameter 
dependence to the susceptibility XpW- I^ R-^f. [lg|, the 
authors succeed in telling apart the classical and quan- 
tum contribution to the Bures metric with consideration 
of these two kinds of dependence. 

The differential of this unitary operator can be divided 
into two parts 

dxU{X)^Di{X)+D2{X), 
^i(A) = Y.^dxu,)P,, D2{X) = Y,u^{^^P^), (14) 



where 9a -Pi = \'<l'i){d\il'i\ + \d\il^i){ipi\- Next, we assume 
that the density matrix p can be simultaneously diago- 
nalized with U , so it can be written as the form 



p = ^PuPi- 



(15) 



The assumption [H, p] = is motivated by the important 
example where both time evolution unitary operators and 
density operators considered are generated by the same 
Hamiltonian i.e., p is the Gibbs thermal state associated 
with H. 
Since 

{PjdxP^)p = 6,jpjji{^,\dx^j,) + {dxiP,\iP,)) = 0, (16) 

after substituting Eq. p4|) into Eq. ^, we obtain 

XpW=4'H^)+x'^'H^), (17) 



xi'HX)^{DlD,),-\{U^D,) 



p\ ' 



(18) 



Here we separate the contributions of Di and I?2- To 
make them explicit, we consider the time evolution oper- 
ator U{X,t) — e"**-^(-^), and assume that the eigenstates 
changing with A can be connected through a smooth uni- 
tary transformation 5(A) which is time-independent and 
satisfies S{X)H{X)S{Xy — i?d(A), where Hd is a diagonal 
Hamiltonian in the fixed A-independent basis. Therefore, 
the unitary operators considered can be written 



u{x,t)^sHx)Ud{x,t)s{x), 



(19) 



with Ud = exp(—itHd). After differentiating with respect 
to A, we get 

Di{X,t) = -itS^{X)[dxHd{X)]S{X)UiX,t), 
D2{X,t)^[A{X),U{X,t)l (20) 

where 

^(A) = idxS{X))^S{X). (21) 

Substituting Eq. ^ into Eq. HH]) leads to 

xi'H^t) ^mdxHdr),- {d,Hd)l]. (22) 



This can be considered as the fluctuation of the quantity 
d\Hd{X) under the state p, with an extra factor t^ . 
The second term of Eq. p?)) will be 

xf(A,t) = ([[/t(A,i),^t(A)][^(A),C/(A,t)])^. (23) 

Since U and p can be diagonalized simultaneously, the 
above expression can be further rewritten as Xp i^it) — 
'^Y.mn Pnn\Amn\'^ (1 " cos[(£'„ - £',„)t]) . This form ex- 

plicitly shows that the time-dependent terms in Xp {^^ t) 
are circular functions. 

Now we would like to make a few general comments 
about Eqs. ^F^ and 



i) They correspond one to one to the two terms in Eq. 
(J12p [with suitably inserted /9„.„]. In particular this 
remark shows that An,m = {En — Ern)~'^ {n\H'\m) . 
Notice that A is nothing but the infinitesimal gen- 
erator of the adiabatic intertwiner mentioned in Sec 

mil 

ii) Xp (i = 1,2) separate apart the contributions of 
the eigenvalue variations from those of the eigen- 
states. In Ref. [16], a similar distinction was made 
for the Bures metric on thermal state manifolds 
(the corresponding terms there were named the 
classical and the quantum, respectively). 



iii) Eqs (P^ and (|^5|) have quite different forms of time 
dependence, xp gives t^ contributions explicitly, 
while the second term consists of circular functions 
with finite period of criticality. When t is large, 
if Xp is not vanishing, it dominates the OFS and 
therefore the decay behavior of the operator fidelity. 

iv) If p is pure [implying, since we have assumed the 
form p = ^,. piiPi, that the initial state is an eigen- 

state of the Hamiltonian] Xp vanishes. Moreover 
if .4(A) commutes with U{X,t), and therefore the 
eigenvectors of U{X) are A-independent, it is obvi- 
ous from ((23)) that Xp vanishes. 



IV. THE TRANSVERSE FIELD ISING MODEL 

In this section, we apply the general formalism devel- 
oped so far to the concrete case of the transverse field 
Ising mode. Before doing that it is useful to discuss a bit 
generally models having a factorization structure. 

A. Factorizable models 

In order to obtain the exact solution of Xp(A,i), we 
analyze the cases where the unitary operator U{X, t) and 
thermal state p(A, /3) have the same factorization struc- 
ture. This assumption hold true, for example, in cases 
where they are both generated by the same Hamiltonian 
H{X). These two quantities can be then written in the 
composite space 

M 

U{X) = (g)t/fe(A), 

fe=0 
M 

pW = (8)Pfe(A), (24) 

fc=0 

where Uk is still an unitary operator and pk is still a den- 
sity operator, corresponding to the fc-th subspace. Note 
that 



M l-l 



M 



dxU{X)=J2[(^UkiX)]®dxUi{X)®[(^ UkiX)]. (25) 



1=0 k=0 



k=l+l 



Since Tr{A (g) B) = Tr(y4)Tr(B) and Tr{pi) = 1, the first 
and second terms of Eq. ([7]) can be expressed respectively 
as 

TT[p{dxU)\dxU)] =Y.TTi[pi{dxUi)HdxUi)] + 5] A,p 

I TrlpUHdxU)] \' = J2\ T^iiPiUHdxUi)] \' +E^«' 

(26) 

wfiere Aw = Tn[piU}{dxUi)\In'[pvUl{dxUv)], and the 
subscript I means the l-th. subspace. 

Thus, the OFS Xp for the factorized unitary operator 

is 

Xp(A) = Exp.'(^) (27) 

I 

Xpj{\) = TYi[pi{dxUi)\dxUi)] 

-\Tn[piU}{dxUi)]\^ (28) 



B. Ising Hamiltoninan 

We now consider transverse field Ising model. The 
Hamiltonian is given by 



M 



H{\) = -J 



l=-M 



(af af+i + Afjf ), 



(29) 



where J is an exchange constant hereafter assumed to 
be the unity, and A is the transverse field strength. This 
model can be calculated analytically by using the Jordan- 
Wigner transformation [22] 



j<i 



(30) 



which maps the spins to fermions. After the Fourier 
transformations, the Hamiltonian in the momentum 
space is 



M 



H{X) 



,2'Kk 



^ [cos( — )-A](44+dV_,-l) 



k=-M 



.27rfc. 



*sin(^)(4dLfe-d_fcdfc), 



(31) 



where N — 2M + 1. This Hamiltonian can be exactly 
solved by a Bogoliubov transformation [23|. However, 
introducing a set of pseudo-spin operators is more conve- 
nient here. Since Uk — n^k {nk — dkdk) commutes with 
every term of the Hamiltonian (|3ip . the pseudo Pauli 
operators can be defined by ^24,] 



S/ca: 



7tjt 



^ky = -^idkdik ~ d^kdk), 
^kz = d\dk + d^_j,d^k - 1- 



(32) 



These give the Pauli matrix in the n^ — n^k = sub- 
space, and become zero matrix in the n^ — ri-fc = ±1 
subspaces. Therefore <^ox, "Joj/i ^Oz are just the standard 
Pauli matrices a^x^ croyj and croz- 

With these operators, the Hamiltonian can be written 
as 

j\/ 

if(A) = ^5l(A)ilM(^)5fc(A) + (A - l)<jo., 
fc=i 

Hk.d — ^k^kz, 



Sk{\) = cxp(-«-— (^fej,). 



(33) 



where 



Vlk = -2\J[\ - cos{2iTk/N)Y + sin2(27rfc/iV), 
2 sin(27rfc/iV) 



Tk = arcsm[ 



^k 



(34) 



We consider the fidelity of two time evolution opera- 
tors U{\) = exp{-itH{X)) and U{X + SX) with ther- 
mal state p — exp{—f3H{X))/Z{P,X), where Z(/3, A) — 
Tr[exp(— /3i?(A))] is the partition function, and f3 — 1/T 
is the inverse temperature. The operators U and p can 
be written in the factorized form 

AI 

U = e-'*(^-i)--®[(g)C/fe], 
fc=i 

p/3(A-l)o-o. Af 



P 



Zo 



[(g) f 



k=l 



Zk 



(35) 



where 



Pk = 5t(A)e-''"'=*^5fe(A) = 6-^*"'=^'^", 
^fen = <;ky sin 6k + <;kz cos 9k , (36) 

and the partition functions are 

Zo = 2cosh(/3(A-1)), Zfe = 2[l + cosh(/3r2fe)]. (37) 



For fc > 0, after substituting Hk^d = ^k'ikz into Eq. ((22|) . 
we have 



xS(A,i)=4t= 



cos" Ok 



cos(/3fife) + f 



(38) 



To calculate X ^j ^^ should first calculate the related 
quantities Ak defined by Eq. (|2ip and the commutator 

[Ak,Uk] = i0 k sin{m k){'ikz sin 6 k -(^ky cos 6k), (39) 

where 0^ denotes dx6k for convenience. Substituting 
them into Eq. (|23p . we obtain 



(2) _ cosh^prik) sin^6>fcsin^(tr?fc) ,.„. 

^p^^ ~ cosh(/3f7fe) + 1 nl • ^ ' 



For fc = subspace, the contribution is 
Thus, we get the OFS 



(41) 



M 



Xp 



(2) 



*=[l-.,„h=(,„A-l,)K4i:.= ^j;g^. 



M 



4E 



2/1 „• 2/ 



cosh(/?rifc) sin 6k sin (tf2 



^ cosh(/?rjfe) + 1 f]^ 



(42) 



Notice that these expressions, for the XY model, can be 
also obtained directly from Eq. (fT2|) . For a given A, we 
can consider the time average OFS which is given by 



1 



XpW = lim 77; / dtXpiX,t). 

1 — >OQ 1 



(43) 




FIG. 2; Time average of the second term of OFS as a function 
of A for a given /3 = 1 with different M. 



Obviously Xp W diverges, since Xp W is proportional 
to t^ . The time average of the circular function contri- 
butions Xp is 



M 



xf\x) = 2Y^ 



k=l 



cos\i{j3Vlk) s\v? Ok 
cosh{l3VLk) + 1 VLl ' 



(44) 



After rescaling Xp ~^ Xp/^j 2TTk/N — > k and taking the 
thermodynamical limit, we get 



9f2 /-^ 
X«(A,t) = — / dk 



X?H\t) = 



1 



■K Jq cosh{(3flk) + 1 



cos^Ok, 



dk 



cosh(/3f2fc) sin^6'fcsin^(irifc) 



IT Jq cosh(/3r2fe) + 1 



^l 



xf\x)=' 



dk 



cosh(/9fife) sin^6'fe 



TT 7o cosh{l3Q.k) + 1 ^l 



(45) 
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FIG. 1: First term of OFS divided by t^ as a function of A 
with different temperatures in the thermodynamical limit. 

One finds that both the quantities Xp ('^j t) ^^^ 
Xp (A) have nontrivial behavior at the critical point 
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FIG. 3: Time average of the second term of OFS as a function 
of A and /3 for a given M = 1000. 



A = 1 and, therefore, can be used as indicators of criti- 
cality. At A = 1 , xp ha-s a minimum (see Fig. [ij , while 
Xp has a maximum, (see Fig. [3]). 

Notice that Xp vanishes only when the inverse tem- 
perature (3 — 00. When the temperature is not zero, 
for long enough time, Xp will be dominated by Xp ^^'^ 
therefore will have a minimum at A = 1. This has to 
be contrasted with the previous results for ground states 
in Ref. ^5] where is was argued that the short time LE 
decay is enhanced at criticality (large N.) For a fixed 

t and A, as temperature increases, Xp grows while Xp 
decays. 

The time average of second term of OFS will diverge 
at the critical points A = 1 under thermodynamical limit 
(see Fig. [5]). It is caused by the denominator flk in 
Eq. (|45p . since it will vanish for those very small ks at 



the critical point A = 1 when N goes to infinity. Thus, 
this divergence will be retained for all /3 (see Fig. [3]). 

In principle, one obtains Xp (A,i) and Xp W from 
experiments in different ways. If we use xp (-^j i) to in- 
vestigate the quantum critical point, the measurement of 
the time interval t is important. However, when xp (-^j t) 
vanishes, we can use Xp W to investigate quantum crit- 

(2) 

icality. This means we can get a Xp (A, t) for a random 
t again and again , then evaluate the time average value. 
By this scheme, we can avoid the demand of the exact 
measurement of the time interval t from the coupling in- 
stant to the measuring instant. 



physical origin and temporal behavior. For the transverse 
field Ising model, we obtained an exact expression for Xp 
and showed that it has nontrivial behavior at the critical 
point both at zero and non-zero temperature. Moreover, 
from p^ it is clear that the OFS depends crucially in the 
level spacing distribution of H. This leads to conjecture 
that OFS might be an effective tool to investigate the 
transition to quantum chaos as well. Finally, we believe 
this type of analysis is directUy relevant to experiments 
aimed at using quantum probes to detect QPTs. 
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